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We construct for the first time, the sequences of stable neutron star (NS) models capable 
of explaining simultaneously, the glitch healing parameters, Q, of both the pulsars, the 
Crab (Q > 0.7) and the Vela (Q < 0.2), on the basis of starquake mechanism of glitch 
generation. These models yield an upper bound on surface redshift of NSs, zji ~ 0.77. 
If the lower limit of the observational constraint of (i) Q > 0.7 for the Crab pulsar 
and (ii) the recent value of the moment of inertia for the Crab pulsar (evaluated on 
the basis of time-dependent acceleration model of the Crab Nebula) , /crab 45 ^ 1-93 
(where I4S = //lO''^ g.cm^), both are imposed together on these models, the models 
yield the value of matching density, Ef, = 9.584 X 10^*gcm~^ at the core-envelope 
boundary. This value of matching density yields a model-independent upper bound on 
neutron star masses, Mmax < 2.22Mq, and the strong lower bounds on surface redshift 
zn ~ 0.6232 and mass M ~ 2.IIM0 for the Crab (Q ~ 0.7) and the strong upper 
bound on surface redshift zn ~ 0.2016, mass M ~ O.982M0 and the moment of inertia 
^Vela,45 — 0.587 for the Vela (Q ~ 0.2) pulsar. However, for the observational constraint 
of the 'central' weighted mean value Q ^ 0.72, and /crab 45 > 1.93, for the Crab pulsar, 
the minimum surface redshift and mass of the Crab pulsar are slightly increased to the 
values zn ~ 0.655 and M ~ 2.149M0 respectively, whereas corresponding to the 'central' 
weighted mean value Q ^ 0.12 for the Vela pulsar, the maximum surface redshift, mass 
and the moment of inertia for the Vela pulsar are slightly decreased to the values zn ~ 
0.1645, M ~ O.828M0 and /yela 45 — 0.459 respectively. These results set an upper and 
lower bound on the energy of a gravitationally redshifted electron-positron annihilation 
line in the range of about 0.309 - 0.315 MeV from the Crab and in the range of about 
0.425 - 0.439 MeV from the Vela pulsar. 
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1. Introduction 

The data on the gUtch heaUng parameter, Q, provide the best tool for testing 
the starquake (Ruderman 1972; Alpar et al 1996) and Vortex unpinning (Alpar 
et al 1993) models of glitch generation in pulsars. Both of these mechanisms of 
glitch generation, in fact, consider NSs, in general, a two component structure: a 
superfluid interior core surrounded by a rigid crust (in the present study we shall 
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use the term 'envelope' which includes the solid crust and other interior portion of 
the star right up to the superfluid core) . In the starquake model, Q is defined as the 
fractional moment of inertia, i.e. the ratio of the moment of inertia of the superfluid 
core, /core, to the moment of inertia of the entire configuration, /total, as (Pines et 
al 1974) 

Q = ^. (1) 

total 

Recently, Crawford & Demiahski (2003) have collected the all measured values of 
the glitch healing parameter Q for Crab and Vela pulsars available in the literature 
and found that for 21 measured values of Q for Crab glitches, a weighted mean of the 
values yields Q = 0.72 ± 0.05, and the range of Q > 0.7 encompasses the observed 
distribution for the Crab pulsar. In order to test the starquake model for the Crab 
pulsar, they have computed Q (as given by Eq.(l)) values for seven representative 
equations of state (EOSs) of dense nuclear matter, covering a range of neutron 
star masses. Their study shows that the much larger values of Q{> 0.7) for the 
Crab pulsar is fulfilled by all the six EOSs (out of seven considered in the study) 
corresponding to a 'realistic' neutron star mass range 1.4 ± O.2M0. By contrast, a 
weighted mean value of the 1 1 measurements for Vela yields a much smaller value 
of Q{= 0.12 ± 0.07) and the all estimates for Vela agree with the likely range of 
Q < 0.2. Thus, their results are found to be consistent with the starquake model 
predictions for the Crab pulsar. They have also concluded that the much smaller 
values of Q < 0.2 for the Vela pulsar are inconsistent with the starquake model 
predictions, since the implied Vela mass based upon their models corresponds to a 
value < 0.5Mq for Q < 0.2, which is too low as compared to the 'realistic' NS mass 
range. 

Thus, in the literature, the starquake is considered as a viable mechanism for 
glitch generation in the Crab and the vortex unpinning, the another mechanism, 
is considered suitable for the Vela pulsar, since it can avoid some other problems 
associated with the starquake explanation of the Vela glitches (see, e.g. Crawford & 
Demianski (2003); and references therein). However, it seems surprising that if the 
internal structure of NSs are described by the same two component conventional 
models (as mentioned above), different kinds of glitch mechanisms are required for 
the explanation of a glitch! Furthermore, it also follows from the above discussion 
that the main reason for not considering the starquake, the feasible mechanism for 
glitch generation in the Vela, lies in the fact that there exists none of the sequence 
of NS models in the literature which could explain simultaneously, on the basis of 
starquake model, both the extreme limiting cases of glitch healing parameter, Q, 
corresponding to the Vela {Q < 0.2) and the Crab {Q > 0.7) pulsars in the range 
< Q < 1.0 for the 'realistic' NS mass values for both the pulsars. The present 
study, therefore, deals with the construction of such models B 



^'however, the other problems associated with the starquake explanation of the Vela glitches (see, 
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We assume that all the NSs belong to the same family of NS sequence which 
terminates at the maximum value of mass. Certainly, this maxima should correspond 
to an upper bound on NS masses. In order to constriict such a sequence, wc have 
to set the extreme causal EOS (in geometrized units), dP/dE = 1 (where P is the 
pressure and E the energy-density) to describe the core. Firstly, because various 
observational studies like - the gamma-ray burst data, X-ray burst data and the 
glitch data etc., and their explanation (see, e.g. Lindblom 1984; Cottam et al 2002; 
Datta & Alpar 1993) favour the stiffest EOSs. The latest estimate of the moment 
of inertia for the Crab pulsar (based upon the 'newest' observational data on the 
Crab nebula mass) rules out most of the existing EOSs of the dense nuclear matter, 
leaving only the stiffest ones (Bejger & Haensel 2002; Haensel et al 2006). Secondly, 
because of the fact that the 'real' EOS of the dense nuclear matter beyond the 
density range ^ lO^^gcm"^ are largely imknown due to the lack of knowledge 
of nuclear interactions (see, e.g. Dolan 1992; and references therein; Haensel et al 
2006), and the various EOSs available in the literature (see, e.g. Arnett & Bowers 
1977) for NS matter represent only an extrapolation of the results far beyond this 
density range. Though, the status of the 'real' EOS for NS matter is not certain, 
one could impose some well-known physical principle, independent of the EOS, such 
as the 'causality condition' {dP/dE = 1) throughout the core of the star beyond a 
fiduciary density, E},, at the core-envelope boundary to ascertain a definite upper 
bound on NS masses (see, e.g., Rhoades & Ruffini 1974; Hartle 1978; Lindblom 
1984; Friedman & Ipser 1987; Kalogera & Baym 1996). In this connection this is 
also to be pointed out here that the maximum mass for any EOS describing the 
core, beyond the density E^, with a subluminal sound velocity turns out to be less 
than that of the upper bound obtained by using the extreme causal EOS (see, e.g., 
Haensel et al 2006). 

The envelope of our models (below the density E}, at the core-envelope boundary) 

may be characterized by the well-known EOS of classical polytrope dlnP/dlnp = Fi 
(where p denotes the density of the rest-mass and Fi is a constant known as the 
adiabatic index) for different values of the constant Fi = (4/3), (5/3) and 2 respec- 
tively. The reason for considering the polytropic EOS for the entire envelope lies 
in the fact that with this EOS, our models yield an upper boimd on NS masses 
independent of the value of Fi, and this upper bound (for a fiduciary choice of i?{,) 
is found fully consistent with those of the values cited in the literature (Kalogera & 
Baym 1996; Friedman & Ipser 1987). Thus, the choice of the said polytropic EOS 
for the entire envelope may be regarded entirely equivalent to the choice of the 
various EOSs like WFF (Wiringa, Fiks & Fabrocini 1988), FPS (Lorenz, Raven- 
hall & Pethick, 1993), NV (Negele & Vautherin 1973), or BPS (Baym, Pethick & 
Sutherland 1971) in an appropriate sequence below the density range E^, adopted 

e.g. Crawford & Demianski (2003); and references therein) are not considered in tiie present paper. 
The future study in this regard may provide some explanation, provided the correlation between 
various parameters of the Crab and the Vela pulsar, obtained in the present study, can be utilized. 
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by various authors in tlie conventional models of NSs (see, e.g., Kalogera & Baym 
1996; Friedman & Ipser 1987), so that the constant Fi appearing in the polytropic 
EOS may be looked upon as an 'average' Fi for the density range below Ei,, speci- 
fied by the sequence of various EOSs in the conventional models of NSs. The choice 
of the constant Fi = 4/3,5/3 and 2 thus become obvious, since this choice can 
cover almost the entire range of density discussed in the literature for NS matter 
which is also applicable for the envelope region - the polytropic EOS with Fi =4/3 
represents the EOS of extreme relativistic degenerate electrons and non-relativistic 
nuclei (Chandrasekhar 1935), Fi = (5/3) represents the well-known EOS of non- 
rclativistic degenerate 'neutron gas' (Oppenlicimmcr & VolkofF 1939), and Fi = 2 
represents the case of extreme relativistic baryons interacting through a vector me- 
son field (Zeldovich 1962) (The value of Fi > 2 is also possible for some EOS 
describing the NS matter, e.g., Malone, Johnson & Bethe 1975; Clark, Heintzmann 
& Grewing 1971, however, the results obtained in this paper remain unaffected for 
the choice of Fi > 2), and the outcome of this study (in terms of explaining the 
glitch healing parameter for various pulsars and predicting the upper bound on the 
compactness of NSs (since the upper bound on mass is independent of the value of 
Fi)) would finally decide, among the chosen values, the 'appropriate' value of Fi 
for the NS envelope. The validity of assuming the extreme causal EOS in the core 
and a polytropic EOS in the envelope of the present models, in view of the various 
modern EOS of dense nuclear matter, is also discussed in the last section of the 
present paper. 

We have noted that in all conservative models of NSs, the choice of the core- 
envelope boundary, rj, (corresponding to a density denoted by Ei,), is somewhat 
arbitrary in the sense that there are no criteria available for the choice of a particu- 
lar matching density, Eh, below which the EOS of the NS matter is assumed to be 
known and unique. One can freely choose somewhat lower values of Ei, (which will 
increase the core size) to obtain higher values of Q (sec, e.g. Shapiro & Teukolsky 
1983; Datta & Alpar 1993). To avoid such a procedure, we choose the core-envelope 
boundary of our models on the basis of the 'compatibility criterion' which asserts 
that for an assigned value of the ratio (a) of central pressure, Pq, to central energy- 
density, Eq, the compactness parameter u{= M/R; total mass to radius ratio) of 
any regular configuration should not exceed the compactness parameter Uh of the 
homogeneous density sphere, in order to assure the compatibility with the hydro- 
static equilibrium (Negi & Durgapal 2001; Negi 2004a). This criterion is capable 
of constraining the core-envelope boundary of any physically realistic NS model. 
A combination of this criterion with those of the observational data on the glitch 
healing parameter and the recently estimated minimum value of the moment of 
inertia for the Crab pulsar (based on the newly estimated 'central value' of the 
Crab nebula mass M(nebula) ~ 4.6M0 in the time-dependent acceleration model), 
-fcrab,45 = 1-93; where /crab,45 = -fcrab/10^® gcm^ (Bejger & Haensel 2003) can 
provide the desired NS models discussed above, since both the theory (criterion) 
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and the observations (stated above) are being used to construct the NS models. 

2. Methodology 

The metric for spherically symmetric and static configurations can be written in 
the following form 

ds'^ = e"dt^ - e^dr'^ - r^dO^ - r^sin^^^^ (2) 

where v and A are functions of r alone. The Oppenheimer-Volkoff (0-V) equations 
(Oppenheimer & Volkoff 1939), resulting from the Einstein field equations for sys- 
tems with isotropic pressure P and energy-density E can be written as 

P' = -{P + E)[4:nPr^ +m]/r{r -2m) (3) 
u'/2 = -P'/{P + E) (4) 
m'(r) = AnEr^ ; (5) 

where the prime denotes the radial derivative and m(r) represents the mass con- 
tained within the radius r 

m{r) = / AivEr^dr. 

Jo 

The coupled Eqs.(3 - 5), are solved for the model (supplemented by the boundary 
conditions: P = E = , m{r = R) = M, e" = e'^ = {1 - 2M/R) = (1 - 2w) at r = 
E) by considering the EOS, dP/dE = 1, in the core and choosing various values of 
Fi in the poly tropic envelope for various assigned values of a such that for each value 
of fj, the compactness ratio of the whole configuration always turns out to be less 
than or equal to the compactness ratio of the corresponding sphere (with the same 
cr) of the homogeneous density distribution and should not exceed the exact absolute 
upper bound on compactness ratio of NSs compatible with causality and pulsational 
stability (see, e.g. Negi 2004b). We find that this condition is uniquely fulfilled by 
all models corresponding to an envelope with Fi = (4/3), (5/3) and 2 respectively, 
if the minimum value of the ratio of pressure to energy-density, Pi,/Eb, at the core- 
envelope boundary reaches about 4.694 x 10^^. The results of the study are presented 
in Fig.l for the value of matching density, Eh = 9.584 x lO^^gcm"^, i.e. about 3.55 
times the nuclear saturation density (note that the particular choice of Eh used 
here turns out to be a consequence of the constraints (i) and (ii) mentioned in the 
abstract of the present study, and therefore, does not represent a fiduciary quantity 
as discussed in the next section). It is seen that the models become pulsationally 
stable up to the maximum value of mass A'/max — 2.22AfQ and radius, R ~ 9.64 — 
10.81 km. The minimum radius results for the model with a Fi = 2 envelope thus 
maximizes the compactness ratio for the stable configuration, u ~ 0.34, as shown 
in Table 1. This behaviour indicates that among various values of Fi chosen in the 
envelope for the density range right from 3.55 times the nuclear saturation density 
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up to zero at the surfaee, the 'average' value of Fi is appropriately described by 
a polytropic index n = 1. This fact will also follow from some special features of 
the model with Fi = 2 envelope discussed in the next section. The binding energy 
per unit mass a[= {Mr — M)/M; where Mr is the rest-mass (see, e.g. Shapiro & 
Tcukolsky 1983)] also approaches a maximum for about 0.3201 for the maximum 
value of mass up to which the configurations remain pulsationally stable. However, 
for the Fi = (4/3) envelope model the binding energy reaches a maximum beyond 
the maximum value of mass. 

3. An application of the models to the Crab and Vela pulsars 

For slowly rotating configurations like the Crab and Vela pulsars (rotation velocity 
about 188 and 70 rad sec~^ respectively) the moment of inertia may be calculated 
in the first order approximation that appears in the form of the Lense-Thirring 
frame dragging-efFect. For Crab and Vela pulsars, the first-order effects turn out 
to be about 1 - 2% (other effects like mass shift and deformation from spherical 
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symmetry due to rotation represent second-order effects which are significant for 
the case of miUisecond pulsars. For Crab and Vela like pulsars, the second-order 

effects turn out to be about 10~^ or even lower; sec. e.g. Arnctt & Bowers 1977; 
Crawford & Demianski 2003. Therefore, these effects can be safely ignored when 
studying the macroscopic parameters of the slowly rotating pulsars as carried out 
in the present paper). These effects are reproduced by an empirical formula that is 
based on the numerical results obtained for thirty theoretical EOSs of dense nuclear 
matter. For NSs, the formula yields (Bejger & Haensel 2002) 

I c^^{l + 5x)MR'^, x>0.1 (6) 

where x is the compactness ratio measured in units of (M0(km)/km), i.e. 

_ M/R _ u 
^ ~ Mg/km ~ 1.477 

only static (non-rotating) parameters of the spherical configuration appear in the 
formula. 

Equation (6) is used, together with coupled Eqs.(3 - 5), to calculate the frac- 
tional moment of inertia given by Eq.(l) and the moment of inertia of the en- 
tire configuration for an assigned value of the matching density, Ei,. It is seen 
that among all models corresponding to different values of Fi in the envelope, 
the lower limits of both the observational constraints for the Crab pulsar, Q > 0.7 
and /crab,45 = -ftotai,45 > 1.93, are appropriately satisfied by the Fi = 2 enve- 
lope model. Since this particular model yields a unique value of matching density, 
Ejy = 9.584 X 10^** g cm~^, for the last two constraints in such a manner that the 
constraint /totai,45 > 1-93 is appropriately fulfilled by all stable NS models corre- 
sponding to the Fi = (5/3) and 2 envelopes for another constraint Q > 0.7 as 
shown in Figs. 2 and 3 respectively (the model with a Fi = (4/3) envelope is ruled 
out since it always yields stable models with Q < 0.7). Another special feature of 
the Fi = 2 envelope model lies in the fact that it yields somewhat wider range of 
glitch healing parameter, < Q < 0.78, for the stable models as compared to the 
Fi = (5/3) envelope model which yields this parameter in the range < Q < 0.75 
for stable models as shown in Fig. 2. Furthermore, the Fi = 2 envelope model yields 
somewhat higher value of maximum surface redshift, zr ~ 0.2016, for the Vela pul- 
sar {Q ~ 0.2) and somewhat lower value of minimum surface redshift, ~ 0.6232, 
for the Crab pulsar {Q ~ 0.7) simultaneously, as compared to the Fi = (5/3) enve- 
lope model. Thus provides somewhat wider range of lower and upper bounds on the 
energy of a gravitationally rcdshiftcd radiation from the surface of the Vela and the 
Crab pulsars respectively. And above all, the compactness maximizes, u ~ 0.34, for 
the Fi = 2 envelope model in the stable sequence of NS models as shown in Fig.l. 
This special feature of the model is also discussed in section 2. 

The above-mentioned features indicate the appropriateness of the Fi = 2 enve- 
lope model and thus emerges an important consequence of this study which states 



(7) 
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Fig. 1. Mass-Radius diagram of the models as discussed in the text for the value of matching 
density E = Ei, = 9.584 X lO"'* g cm^'^ at the corc-cnvclope boundary. The labels a, b and 
c represent the models for an envelope with Fi = (4/3), (5/3) and 2 respectively. The minimum 
value of the ratio of pressure to energy-density, (Pj,/i5(,), at the core envelope boundary is obtained 
as 4.694 x 10^'^, such that for an assigned value of a, the inequality u < u/j is always satisfied for 
all the models as shown in Table 1. 



that the envelope of 'real' NSs may be well approximated by an EOS of classical 
polytrope with 'average' value of polytropic index, n, closer to 1. 

For the minimum value of Q ~ 0.7, Fig. 3 yields /total, 45 = /crab.45 — 1-93 for 
the Fi = 2 envelope model and by the use of Fig. 2 and Table 1, wc obtain the 
minimum values of mass, M ~ 2.IIM0, and surface redshift, zr ~ 0.6232, for the 
Crab pulsar. On the other hand, for the maximum value of Q ~ 0.2 which belongs to 
the Vela pulsar, Fig. 3 yields /total. 45 — 0.587 corresponding to the Fi = 2 envelope 
model and Fig. 2 and Table 1 yield the maximum values of mass, M ~ O.982M0, 
and surface redshift, zr ~ 0.2016, for the Vela pulsar (note that the upper weighted 
mean value of Q « 0.19 for the Vela pulsar. Fig. 2 and Table 1 (for the Fi = 2 
envelope model) yield the maximum value of mass, M ~ O.96IM0, and surface 
redshift, zr ~ 0.1966, for the Vela pulsar. Obviously, these values are much closer 
to the values corresponding to the maximum value of Q » 0.2 for the Vela pulsar). 
However, for the observational constraint of the 'central' weighted mean values 
of Q ^ 0.72 for the Crab and Q « 0.12 for the Vela pulsar, Fig.2 and Table 
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Fig. 2. Fractional moment of inertia Q{= /corc//total) vs. total mass M for the configurations 
presented in Fig.l. The labels a, b and c represent the models for an envelope with Fi = (4/3), (5/3) 
and 2 respectively. 

1 (for the Fi = 2 envefope model) yield the slightly increased values of mass, 
M ~ 2.149M0, and the surface redshift, zn ~ 0.655, for the Crab and somewhat 
decreased values of mass, M ~ O.828M0, and surface redshift, zj^ ~ 0.1645, for the 
Vela pulsar respectively. For these central weighted mean values of Q, Fig. 3 yields 
for the Fi = 2 envelope model, the moment of inertia /crab,45 — 1.971 for the Crab 
and /vcia,45 — 0.459 for the Vela pulsar respectively. For the lower weighted mean 
value of Q w 0.05 corresponding to the Vela pulsar. Fig. 2 and Table 1 yield (for the 
Fi = 2 envelope model) the minimum values of mass, M ~ O.685M0, and surface 
redshift, zr ~ 0.1312, whereas Fig. 3 yields a minimum value of moment of inertia, 
A^eia,45 — 0.355, for the Vela pulsar. 



4. Implications of the models for the extraordinEtry gemima-ray 
burst GRB790305b 

Apart from the models of the Crab and Vela pulsars, let us consider two of the main 
findings related to the extraordinary gamma-ray burst of 5 March 1979: (i) it gives 
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Fig. 3. Fractional moment of inertia Q{= /corc//total) vs. moment of inertia of the whole strueture 
/total, 45 for the configurations presented in Fig.l. The labels a, b and c represent the models for 
an envelope with Fi = (4/3), (5/3) and 2 respectively, /total, 46 is defined as /total / 10"*^ gcm^. 



the only reliable estimate of the surface redshift, zr = 0.23 ±0.07 ( after taking due 
account of thermal blueshift), associated with the supernova remnant N49 in the 
Large Magellanic Cloud (sec, e.g. Higdon & Lingcnfcldcr 1990: Douchin & Haensel 
2001; and references therein), and (ii) the implied peak luminosities of the repeating 
burst GRB790305b correspond to an energy of 10^* ergs, which is possible only when 
a starquakc releases at least 10~® of the NS gravitational binding-energy of 10^^ 
ergs (Higdon & Lingenfelder 1990). In order to reproduce both of these findings, one 
would require a NS model compatible with starquake model predictions that could 
also account for a surface redshift zr = 0.23 ± 0.07. Both of these requirements 
are, in fact, fulfilled by our models. For the case of the Fi = 2 envelope model, 
we get from Fig. 2 - 3 the central value of surface redshift zr = 0.23 for a mass 
M ~ 1.O95M0 with Q ~ 0.259. The binding energy corresponding to this case is 
obtained as 2.618 x 10^^ ergs which is capable of releasing 2.618 x 10*^ ergs of energy 
required for the latter burst. 
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5. Results and conclusions 

This study constructs the stable sequences of NS models terminate at the value of 
maximum mass, Mmax — 2.22M0, independent of the EOSs of the envelope, for 
the matching density, Eb = 9.584 x 10""^^ g cm~^, at the corc-cnvolope boundary. 
This value of 'matching density' is a consequence of the observational constraints 
Q ~ 0.7 and /crab,45 — 1-93 (associated with the Crab pulsar) imposed together on 
the Fi = 2 envelope model and in this sense does not represent a fiduciary quantity. 
The upper bound of the surface redshift, zn ~ 0.77 (corresponding to a u value 
~ 0.34), however, belongs to the model with a Fi = 2 envelope which is consistent 
with the absolute upper bound on the surface redshift of NS models compatible 
with causality and pulsational stability (Ncgi 2004b). This special feature, together 
with some other remarkable ones, discussed in the present study underline the 
appropriateness of the Fi = 2 envelope model. 

Since among the variety of modern EOSs discussed in the literature, the upper 
bound on NS mass compatible with causality and dynamical stability can reach 
a value up to 2.2M0 (in this category, the SLy (Douchin & Haensel 2001) EOS 
yields a maximum mass of 2.O5M0, whereas the BGNl (Balberg & Gal 1997) and 
the APR (Aknial ct. al. 1998) EOSs yield the maximum masses of 2.18Af0 and 
2.2IM0 respectively (see, e.g. Haensel et al 2006)). In view of this result, the model- 
independent maximum mass, M^ax — 2.22M0, obtained in this study may be 
regarded as good as those obtained on the basis of modern nuclear theory. 

In addition to this result, the appropriate sequences of stable NS models obtained 
in this study can explain the glitch healing parameter, Q, of any glitching pulsar, 
provided the weighted mean values of Q lie in the range < Q < 0.78. This finding 
also reveals that if the starquake is considered to be a viable mechanism for glitch 
generation in all pulsars, then the envelope of 'real' NSs may be well approximated 
by a polytropic EOS corresponding to a polytropic index, n, closer to 1. 

For the value of matching density, Ef, = 9.584 x 10^* g cm~^, the Fi = 2 envelope 
model yields the minimum values of mass M ~ 2.llMr.) and surface redshift z/? ~ 
0.6232 for the Crab {Q ~ 0.7) and the maximum values of mass M ~ O.982M0 
and surface redshift zn ^ 0.2016 for the Vela pulsar (Q ~ 0.2). The minimum 
mass and surface redshift for the Crab pulsar are slightly increased up to the values 
M ~ 2.I49M0 and z^ — 0.655 respectively, if the 'central' weighted mean value 
of <5 ~ 0.72 and the moment of inertia /crab,45 > 1-93 are also imposed on these 
models. However, for the 'central' weighted mean value of Q ^ 0.12 corresponding 
to the Vela pulsar, the maximum mass and surface redshift are somewhat decreased 
to the values M ~ O.828M0 and zr ~ 0.1645 respectively. This value of mass 
and surface redshift for the Vela pulsar can further decrease up to the values M ~ 
O.685-M0 and zn, — 0.1312 respectively, if the lower weighted mean value of Q ~ 0.05 
for the Vela pulsar is imposed. These results predict the upper and lower bounds on 
the energy of a gravitationally redshifted electron-positron annihilation line in the 
range of about 0.309 - 0.315 MeV from the Crab and in the range of about 0.425 - 
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0.439 MeV from the Vela pulsar respectively. 

For a comparison, if the observational constraint of the minimum value of 
J'crab.45 — 3.04 (the valuc of moment of inertia for the Crab pulsar obtained earlier 
by Bejger & Haensel (2002), on the basis of the constant- acceleration model for 
the Crab nebula) together with Q ~ 0.7 is imposed on the models studied in the 
present paper, the Fi = 2 envelope model yields the value of matching density. 
El, = 7.0794 X 10^"* g cm"'^. This value of Ei, yields a modcl-indcpcndent upper 
bound on NS mass Mmax — 2.59M0. This value of maximum mass, however, repre- 
sents an 'average' of the maximum NS masses in the range 2.2M0 < Mmax < 2.9M0 
obtained by Kalogera & Baym 1996 (and references therein) on the basis of other 
EOSs for NS matter, fitted to experimental nucleon-nucleon scattering data and 
the properties of light nuclei. For this lower value of matching density, the Fi = 2 
envelope models yield the minimum value of mass M ~ 2.455A'/0 for the Crab 
{Q ~ 0.7) and the maximum value of mass M ~ 1.142M0 for the Vela pulsar 
{Q ~ 0.2). The minimum mass for the Crab pulsar is slightly increased up to the 
value M ~ 2.5M0, if the 'central' weighted mean value of Q ~ 0.72 and the moment 
of inertia /crab.45 > 3.04 are also imposed on these models. However, corresponding 
to the 'central' weighted mean value of Q ~ 0.12, the maximum mass of the Vela 
pulsar is somewhat decreased to the value M ~ O.964M0. This value of mass for 
the Vela pulsar can further decrease up to the value M ~ O.796M0, if the lower 
weighted mean value of Q ~ 0.05 for the Vela pulsar is imposed. 

Furthermore, the study can also explain some special features associated with 
the extraordinary gamma-ray burst of 5 March 1979. 
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